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C^*l . We discuss the inspiral of a small body around a Kerr black hole. When the time scale of the radiation reaction 

is sufficiently longer than its orbital period, the leading order orbital evolution is described only by the knowl- 
edge of the averaged evolution of the constants of motion, i.e., the energy, azimuthal angular momentum and 
the Carter constant. Although there is no conserved current composed of the perturbation field corresponding to 
■ the Carter constant, it has been shown that the averaged rate of change of the Carter constant can be given by a 

D ' simple formula, when there exists a simultaneous turning point of the radial and polar oscillations. However, an 

, inspiralling orbit may cross a "resonance" point, where the frequencies of the radial and polar orbital oscillations 

are in a rational ratio. At the resonance point, one cannot find a simultaneous turning point in general. Hence, 
even for the averaged rate of change of the Carter constant, a direct computation of the"self-force", which is 
still challenging especially in the case of the Kerr background, seems to be necessary. In this paper, we present 
a method of computing the averaged rate of change of the Carter constant in a relatively simple manner at the 
O ' "resonance" point. 

cr 

5^ , PACS numbers: 04.30.Db, 04.25.Nx, 04.20.Cv 

• I. INTRODUCTION 

in 

CO 

f^*) ' The study of the two-body problem in general relativity has been attracting much attention in recent years to obtain a better 
understanding of the gravitational wave emission from a binary system. If one body is much less massive than the other, i.e., 
' when we consider extreme mass ratio inspirals of a compact object into a super massive black hole, the compact object can be 
treated as a point particle and the two-body problem is handled within the self-force program (See reviews Jl] |2[ and references 
£f) • therein.) 

When a nonspinning particle with mass ^ is moving along a bound orbit on a Kerr black hole with mass M(^> fj,) and the 
spin parameter a, the long-term evolution of the orbit due to the radiation reaction of the gravitational wave emission is dictated 
• rH | by the long-time averaged rates of change of three constants of motion: the energy E, the azimuthal angular momentum L and 
the Carter constant Q y[] since the characteristic time scale of the secular orbital evolution is sufficiently longer than the orbital 
period. The orbital phase errors originating from the other effects scales as 0((M/ fj,) ). This particular approximation of the 
orbital evolution is referred to as the adiabatic approximation JJIst]- 

On average, the losses of the energy and the azimuthal angular momentum of the particle are to be balanced with those carried 
by gravitational waves to the infinity or into the Kerr horizon thanks to their global conservation laws. This is nothing but the 
well-established balance argument. Although we do not have any balance argument for the Carter constant, based on Mino's 
argument JH, Sago et al. iflolfTTll (see also Drasco et al. |[l2tl for the scalar case) have succeeded in deriving a formula for the 
averaged rate of change of the Carter constant (dQ / dt) similar to the one obtained by the balance argument for the energy and 
azimuthal angular momentum. 

However, the above mentioned work assumes that the particle's radial and polar orbital frequencies, fi r and fig, are in an 
irrational ratio. When O r and Slg cross a resonance point, at which £l r /flg = j r /je holds with coprime integers j r and jg, the 
above formula is not correct, and the phase error due to the ignorance of the correct formula can be as large as 0((M/ /j,) 1 / 2 ). This 
was first pointed out by Mino (4\ (See also Ref. JH). Recently, Flanagan and Hinderer lfl3ll investigated this issue quantitatively 
within the post-Newtonian approximation. For the inspiral in the extreme mass ratio regime M ^> [i, this phase error is much 
larger than the one from the first order conservative self-forces or the long-time averaged second order dissipative self-forces, 
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which scale as 0((M//x)°) JH-0, H-QjJ . For this reason, it is very demanded to derive a formula for (dQ/dt) in the adiabatic 
approximation for a resonance orbit. 

The main differences from the off -resonance case in deriving (dQ/dt) are in the following two points. First, the orbit is 
no longer characterized by the three constants of motion alone. Instead, as was pointed out by Flanagan et al. 11711 . the orbit 
depends on the relative offset phase specified by AA, the separation between the instances when the r- and 8- oscillations reach 
their minima. The offset phase also slowly evolves compared to the orbital period near resonance points. Thus, the resonance 
orbit must be characterized by four parameters: {E,L Z ,Q, AA}. Second, we cannot use the Mino's argument 10], which proves 
that (dQ /dt) can be evaluated by using the radiative field defined by "the half -retarded minus half-advanced field", which is by 
definition regular at the location of the particle. The key ingredient of the Mino's argument is the invariance of the geodesies 
under the transformation (t, r, 8, (j>) — > (— t, r, 8, —(f)) that is assured if and only if the r- and 8- components of the body's four 
velocity simultaneously vanish at some time, i.e., when we can choose AA = 0. Clearly, this condition is not satisfied for generic 
resonance geodesies. 

Our final goal is to obtain the evolution of a resonant inspiral with the phase accuracy better than unity. Toward this ambitious 
goal, we here begin with deriving a practical formula of (dQ/dt) for resonant inspirals. In the resonance case, we find that 
(dQ /dt) acquires both the contributions from the radiative and symmetric fields, the latter defined by "the half-retarded plus 
half-advanced field". In general, the symmetric field diverges at the particle's location, thus requires some regularization. In 
the later section, we will propose a new regularization scheme for the symmetric portion of (dQ/dt). The key point of our 
regularization method is that it is compatible with the Teukolsky formalism lfl8ll . This makes the calculation much simpler than 
the self-force calculation, which is already numerically challenging for a particle orbiting the Schwarzschild black hole (see 
Ref. (H] for recent progress in this topic). 

In this paper, to avoid complications, we focus on a point scalar particle coupled to a massless field with the mass fi and the 
charge q. Although there seems no essential obstacle in extending our scalar formalism to the gravitational case in principle, the 
latter is more involved since it requires the reconstruction of the metric perturbations, which is a little complicated in the Kerr 
case lfl9l|20tl. 

The remainder of this manuscript is organized as follows. In Sec. |IIJ we begin with describing the model of a point scalar 
particle coupled to a massless scalar field. Then, we review the generic bounded geodesic motion around a Kerr black hole, 
particularly focusing on the resonance case. The retarded solution of the scalar field equation expressed in the form of the 
mode decomposition is also summarized there. Then, we derive a general formula for the long-time averaged rate of change 
of the Carter constant in Sec. [HI] We separate the formula into two pieces, the radiative and the symmetric parts that originate 
from the radiative and symmetric fields, respectively. The result for the radiative part is shown in Sec. IIIIBI and the result for 
the symmetric part is presented in Sec. IIII CI where we also discuss the mode sum regularization of the symmetric part since 
it diverges at the particle's location. Finally, we summarize this paper in Sec. [IV] To clarify the impact of our work, we also 
present there how the long-term evolution of the resonance orbit is described in the adiabatic approximation. In Appendix [Al 
we discuss the possibility that the resonance effect on the phase error mig ht be larger than 0{{n/M)~^ 2 ). Some detailed 
derivations related to the regularized symmetric part are summarized in Appendix IE1 

Throughout this paper, we use geometrical units G = c = 1 and the sign convention of the metric is (—,+,+,+). We 
basically follow the notation used by Drasco et al. lfl2ll . 



II. PRELIMINARIES: THE MOTION OF A POINT SCALAR PARTICLE 

We consider a point scalar particle with the bare rest mass fi = 1 1 and the scalar charge q, moving in the Kerr spacetime. The 
Kerr metric in the Boyer-Lindquist coordinates is given by: 

, „ , „ ( 2Mr\ , 2 4Marsm 2 8 , ,, £ , 2 „ I/l2 f ? ? 2Ma 2 r . 2 A . 2 „ , , 2 
g^dx^dx" = - ( 1 — j dt 2 dtd<j) + —dr 2 + Y,d8 2 + ( r + a H — sin 2 8 J sin 2 8 def, (1) 

where M and a are the mass and the Kerr parameter, respectively, and 

£ := r 2 + a 2 cos 2 8, A := r 2 - 2 Air + a 2 . (2) 
There are two Killing vectors associated with the stationarity and axisymmetry of the Kerr spacetime: 



Strictly speaking, fi is not constant because of the self-force effect. However, the correction due to this time dependence is 0((M/ in the orbital phase, 
i.e. suppressed as well as the conservative self-force correction. Since we focus on the correction of 0((M/ /i) 1 / 2 ), we neglect the change of the mass in this 
paper. 
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The Kerr spacetime also has an irreducible Killing tensor given by 

K^ v := 2EZ( M 7i„) + r 2 g^ = 2Em ((J %) - a 2 cos 2 9g^, (3) 

with the Kinnersley null tetrad defined by P := (r 2 + a 2 , A, 0,a)/A, := (r 2 + a 2 , -A, 0, a)/(2E) and := 
(iasm9, 0, 1, i/ sin6)/(y/2(r + ia cos 6)) where l^n^ = — 1 and m^fh^ = 1 while all the other inner products vanish. 



A. Bounded geodesic motion around a Kerr black hole in resonance 

Neglecting the effect of the self-interaction, a particle can be regarded as a test particle that moves along a geodesic on the 
background spacetime, 2^(t) := (i(r), r(r), 0(t), 4>{t)) parametrized by the the proper time r. There are three integrals of 
motion owing to the symmetries of the Kerr spacetime: the energy E, the azimuthal angular momentum L and the Carter constant 
Q 10], respectively expressed as 



E := -u a & = / ! _ __j „t 



2Mr\ t 2Marsm 2 9 



- II: ■ 



E 

i : = ^ =- 2MQ ; sin2g ^+ (r2+a2) ; Aa2sin2g sin 2 ^^ 

Q := u a u K aB = ( - a sin 9E ) + a 2 cos 2 + E 2 (u e ) 2 , (4) 
\sin# J 

where u a := dz a /dr is the particle's four velocity. (Note that E and L are 0(/i), and Q is 0(/i 2 ) if we recover /j,.) 

For later analysis, it is convenient to parametrize the particle's orbit as z(X) := (t(X), r(X),9(X), <fi(X)) with the "Mino time" 
A, which is related to t as dX := dr/Y; fl. Using A, the equations of motion are written as 

£) 2 = «<->• (^)" 

^ = VJe(0) + Vk(r) , ^ = V^(tf) + VV(r), (6) 



where 



R(r) := P(r) 2 - A(r 2 +Q), 

P(r) := E(r 2 + a 2 ) - aL, 

e(cos6>) := C-{(7W(l-£ 2 ) + L 2 }cos 2 6> + a 2 (l-£ 2 )cos 4 

Vitf(fl) := -a(aEsin 2 9-L), V tr (r) := L-±^-P(r), 



-a£+-4-, ^ r (r) := £p(r), 

sm 6 1 A 



with C := Q — (aE — L) 2 . Here it should be noted that the equations of r and 9 are completely decoupled, and that those of t 
and 4> are divided into r and 9 dependent parts. 

First, it follows from Eqs. (0 that the r- and 9- oscillations are independently periodic for a bound orbit: 

r(A) =r(A + 7iA r ), 9{X) = 9(X + kAg), (7) 

where n and k are integers and the periods A r and Ag with respect to the Mino time A are given by 

with the turning points r m j n , r max and # m i n < 7r/2. For later convenience, we introduce f and 9 which denote the fiducial 
solutions of Eqs. (0 such that r(nA r ) = r m j n and 9(kA r ) = m ; n . 
The precise meaning of the resonance is that the orbital frequencies 

2tt ^ 2tt 
T r := — , T e := — 
A r A 
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are related to each other as 

£ - 5, (9) 
where j r and jg are coprime integers. We here define the resonance frequency, T, and the corresponding period, A, as 

jr je T 

In the resonance case, in general, there is a difference between the Mino times of reaching the minima of r- and 8- oscillations. 
We call this difference "offset phase" and denote it by AA. Note that AA can be shifted by 

A' := — 

irje 

depending on which minima we focus on. Therefore, we choose A A such that |AA| < A'/2. Setting the origin of the A 
coordinate to satisfy 0(A) = 0(A), the orbit is written as 

r(A) = f (A — AA) , 0(A) = 0(A) (10) 

without loss of generality. 

Second, the motion in the t- direction is also separated into r- and ^-dependent parts. We introduce the averaged values of 
V tr (r) and V tg (9) defined by 

1 /-A-r 1 fAs 

-L / iW Tr r_/w\n i-w T \ 



(Vtr)~T dX' V tr [f(X% ( Vte ):=— d\' V t g[6(X% (11) 
A r Jo A e Jo 

In integrating the t-component of the equation of motion, the choice of the origin of the ^-coordinate can be ch ang ed freely 
because of the time-translation symmetry of the Kerr spacetime. Then, without loss of generality, t(X) is given by 11211 

t(X) = YfA + At r (X — AA) + Ate(A), (12) 

where T t := (V tr ) + (V t g), and 

At r (X):= f dX'{V tr [f(X')} - (V tr )}, At g (X):= f dX'{V t g[6(X')} - (V tg )}. (13) 
Jo Jo 

The function 0(A) is deduced in the same way as t(X). Again considering the axial symmetry of the Kerr spacetime, we have 

0(A) = T^A + Acj) r (X - AA) + A(f>s(X) , (14) 
where we define := (V^.) + (V<pg} and 

A0 r (A):=/ dX'{V^ r [f(X')}~ (V^ r )} , A0 e (A):= f d\'{V&[§(\')]- {Vp)} (15) 
Jo Jo 

with 

i rK -, Me 

J- / r_/. /m It -r \ J- 



{V^)'.= -r-\ dX'V^ r [f(X% (V^g):=— dX'VteW)]. (16) 

A r Jo A <? Jo 

From the above discussion, the orbit is expressed as 

z(X) = z(X — AX, A) (17) 

with the two variable function 

z(X r , Xg) = [T t Xg,0, 0, T A e ] + Az(A r , Xg), (18) 
where the oscillating part of the orbit is denoted by 

Az(A r , Xg) = [At r (X r ) + Atg(Xg),f(X r ), 6(X e ), A0 r (A r ) + A0 e (A e )]. (19) 
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B. The retarded solution of the scalar field equation 

The equation of the scalar field induced by a charged particle is given by 

«rV,V„$(i) = 7(x), T(x) := -q [ +0 ° dT S(4 ] ix ~ Z }^\ (20) 

J-oo V-5(z(r)J 

where V M and g show the covariant differentiation and the determinant for the background Kerr metric, g^ v , in Eq. (Q]), respec- 
tively. (See Ref. |[l2tl for further details). Assuming that there is no incoming scalar wave, we take the retarded solution, $( ret ). 
To obtain the solution, we here construct the retarded Green function in terms of the mode functions of the field equation. 

From the fact that the scalar field equation is separable in the Kerr spacetime lfl8[l2ltl . we can write the mode functions in the 
form of 

nl lm (t, r, 9, 4>) := J—^ e~^ ^L»(r*)Wfl, <t>), (21) 
V t + or 

where r* is the tortoise coordinate defined by dr* := ((r 2 + a 2 ) / A)dr, and S tJ i m (9, 4>) := < d U jim{&)e %m ^ is the spheroidal 
harmonics normalized as 

sin 8S u i m {8,4>)* Sui> m > {B,4>) =8 lv 8 mm i. (22) 

The superscript b represents one of the four distinct boundary conditions: "up", "down", "in" and "out". The corresponding 
radial functions are respectively defined by lfl2l I22I 

T^i m \p m u\- 1/2 e- l P- r ' , (r*^-oo), 









— Pwlm 


<Tm(r*) 


- 7/ in * 

— u ujhn 

up * 
ujlm 



-1/2 [ e -*-r- +aulme iur^ ; ( r *^oo), 



|Pmcj I 



-1/2 



, (t* — > —00), 



2 e iwr *, (r*^oo), 



(23) 



with p TOW := w - ma/ (2Mr + ) and r + := M + \JM 2 ~ a 2 . The coefficients a^im and /3 w ; m are normalization constants. T^im 
and a u i m are, respectively, complex transmission and reflection coefficients. Both the complex coefficients /U W ; TO and are 
also defined by Eq. d23l . In this work, we choose the normalizations of the "in" and "up" solutions as 

a ulm = a -io I —mi fiulm = '— w I— mi (24) 

so that the invariance of the radial equation under the transformation (u, m) — > (—a;, — m) implies the identity 

«-*j-m(0 = ulUr*)- (25) 

From the identity d25l l, we find 

T u>lm ~ T -uil-m, a ujlm ~ V—wl—mi H'ulm ~ H*—ul—m-> v u>lm ~ v —ul—m- 

It should be also noted that we choose the phase of the normalization of the angular function a) ; m (6') so that 

9- w J-m(*) = (-l) m e wim (0). (27) 
The retarded Green function for the scalar field equation (l20l is expressed as the factorized form in terms of the mode functions 



f + 00 +oo I 

G(rCt) ^' x ') = 77^7 / duj T, E T7TT „ r KU^L'^Wr r') + 7^)7^ *{x')H{r' r)] 

1D7TZ J_ oc |uJj CKLulrnPLulrn 

1—0 m——l 
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KUr)Kll*(r')H(r r') + RZ? m * (r)R u J lm (r')H(r' r) 



(28) 



with the Heaviside function H(x) :— j^ oo S(y)dy, and R^i m {r) = 2w^ m /\/r 2 + a 2 . In showing the second equality of 

Eq. (1281) . we used the relations between the radial functions u^ lm (r*) in Eq. (1231 . This Green function leads the retarded 
solution as 



^\ x ) = / d i x'yf^)G {Te%) {x,x')'J{x') = -q / drG (rot) [x, z(t)] . 



(29) 



It is helpful for the later analysis to decompose the retarded field into a linear combination of the radiative (anti-symmetric) 
and symmetric fields as 



$ (rct) (x) := $ (rad) (x) + ^ symJ (x), 



(sym) / 



$ (rad) (x) = -q / drG (rad) [x,z(r)], $ (sym) (x) = -q / drG (sym) [x,z(r)] 



where 



with 

G (rad) (x,x') = i [G (rot) (x,x')-G (adv) (x,x') 
G (sym) (x, a;') = i [G (rot) (a, x') + G (adv) (x, a;') 

Here the advanced Green function, G ( - adv ^ (x, x'), is derived from the retarded Green function as 

G (adv) (x,x') = G (rct) (x',x). 
The radiative Green function can be rewritten as lll2"ll 

X3 +OO I 



(30) 
(31) 

(32) 
(33) 

(34) 



G (rad) (*, x') = J- 

327Ti 



E F7 

i=0 m=-J 1 1 



out / \ out*/ l\ , WPmu \7ujlm\ down/ \ down*/ /\ 



(35) 



and the Heaviside functions are eliminated. On the other hand, we rewrite the symmetric Green function in the following way. 
Writing G' adv -'(x, x') in terms of the mode functions in a similar manner to Eq. (ESI , replacing (oj,m) to (— w, — m), and 
rearranging it with Eqs. (l24l . (l25t . ( l26l . and (ITTt . we find 



G (adv) (x,x') 



167ri 



+00 +00 l 



1=0 m=- 1 ^ a ulmPu,lr, 



x [K P MKT m (r')H(r r) + - r) 

Substituting Eqs. d28l and d36t into Eq. d33l l, we find that the symmetric Green function is rewritten as 



, p+00 +°° ( 

g^Hx, x 1 ) = / ^x; e ^rWo, 



(36) 



(37) 



where 



sSTW) = e--(*-*')+*™w-*')e <1 , Im (fl)e: Jm ((9') 



1 



if, 



KUr)KT m *(r')H(r - r>) + (r)R% m (r>)H(/ r) 



(38) 



and represents the imaginary part of X. 
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III. ADIABATIC EVOLUTION OF THE CARTER CONSTANT AT RESONANCE 

We denote the three constants of motion as P := {E, L, Q}. The time averaged rates of change of P can be reexpressed as 

'dP\ 1 I dl 



- 



dt I , T t \ dX i A 



(39) 



with the definition of the A average, 



1 



= lim ^= [ +1 dX'F(X'). (40) 



One key property of the adiabatic approximation developed in Ref. |4J] for off -resonance orbits is that (dP jdX^j can be evaluated 
by using the radiative solution of the field equation ( 1201 as 



where / Q is a differential operator defined by II23I - I25H 

/ Q [$(x)] := g( 5 a ^ + u a u f3 )\/ p ^(x). (42) 

Since the radiative field has no divergent behaviour at the location of the particle, this relation gives (dP /dX) . without any 
regularization procedure. The final expressions for (dP /dX) . are given in terms of the asymptotic amplitudes of the scalar field 
at infinity and on the horizon lfl2ll . 

The derivation of Eq.fiTb relies on the assumption that there exists a simultaneous turning point of the r- and 9- oscillations 
of the particle. In the case of a resonance orbit, however, the frequencies of the r- and 9- oscillations are in a rational ratio, and 
hence one cannot find a simultaneous turning point in general. This means that the orbit depends on the non-vanishing offset 
phase, AA. Therefore, the expressions for (dP/dt) t should be derived without using the relation (|4TV 

As for (dE/dt) t and (dL/dt) v the balance argument can be used to relate them to the fluxes of the energy and angular 
momentum to the infinity or into the black hole horizon. Therefore, one can conclude that the expressions for (dE/dt) t and 
(dL/dt} t in the resonance case are the same as those in the non-resonance case fl7il . On the other hand, there is no known flux 
composed of the perturbation field corresponding to the Carter constant, and hence we cannot use the balance argument to derive 
the averaged rate of change, (dQ/dt) t . We need to reformulate the calculation of (dQ/dt) t without using the relation ( T4"TV 
In this section, we derive a formula for computing (dQ/dt) t in the resonance case, by starting from the basic formula of the 
self-force acting on a particle. 



A. Formal reduction of (dQ/dt) t for a resonance orbit 

In contrast to (dE/dt) t and (dL/dt) t , we need to directly evaluate the time derivative of the Carter constant defined in Eq. (0]i 
caused by the self-force acting on the particle, which is given by 

^ = 2EA>u a a' 3 (43) 

with the aid of the Killing tensor equation V( 7 A" Q(3 ) = 0. Here a a is the four acceleration due to the scalar self-forces with the 
retarded scalar field given by [26] 

a a [z( T )] := lira f a [^ (R \x)}, (44) 

where & R > is the so-called i?-part of the retarded solution, which satisfies the homogeneous equation of Eq. (1201 . is 
schematically expressed as 

$W(x) := $ (rct >(x) ~& s \x) 

= $( rad )(x) + {V sym >(2;) - $ (5) (:c)} , (45) 

where (x) is the S'-part of the retarded field j26ll . which satisfies the same equation as the retarded (or symmetric) solution. 
It should be noted that, the difference between $( s y m ) and $( s ) is regular at the location of the particle while each of them is 
singular. 



Now plugging Eq. d44b into Eq. ( 1431 , the time derivative of the Carter constant is rewritten in terms of the i?-part of the 
retarded field ^^(x): 



dX dX 



(46) 



Since the first term is a total derivative up to the leading order in q, it does not contribute after taking a long-time average. Hence, 
we omit the first term in Eq. (|46V Substituting the explicit expression for K a p into Eq. (l46l i. and dropping the first term, we 
obtain 



dQ\ J2q 
dt /, ' T t 



[6 X [z(x)& R \x)]] 



x—z{\) j \ 



(47) 



where X is a differential operator defined by 

6, := V tr {r)d t + V^(r)fy + ^T^. 
Following the decomposition of the retarded field in Eq. (1451 . we also decompose (dQ/dt) t into two parts as 



dQ\ /^\ (rad) + /^\ 
It \ dt /t \ * A 



(sym— 5) 



with 



dQ\ (rad) 



2g 2 ,. 1 



r/A 



dt j 



(sym-S) 



— T 
T 



d\' \6 X h(x)^(x')G^ ad \x,x')\) _ 



(48) 



(49) 



(50) 



x=z{X),x'=z{X') 



(51) 



where G' 5 - 1 (x, x') is the 5-part of the retarded Green function (see Eq. (|8~H below). 

We reduce Eqs. d50b and ( BTT i into more tractable expressions. For this purpose, we here consider a function of two variables, 
G(x,x') = T,{x)H(x')G{x,x'). Using the identity, 



^-G(z(X r ,Xe),x') 



{-lxV dr + ~^xT dt + dX~ r ' G(X ' X ) 



(52) 



x=z(X r ,\e) 



we find 



Jim — / dA {6* *')!} m 
= y dA e dA r <5(A r -A e + AA){o, [g^z')]^ 



^(A„A e ) 



+ 00 



dA e / dX r 6{X r -X e + AA) «V tr )d t + G(x,x') 



x— 2(A r ,Ae) (iA r 



G (z(A r , Ag), a;') 



((Vir)a t +(^)G(l,l') 



x=z(\) / x d(AX) 



(G[z(X),x']) ; 



(53) 



Here we insert the <5-function 6(X r — Xg + AA) on the second line to deal with the r- and 0-oscillations separately, keeping the 
information of the offset phase. In the last equality, we performe integration by parts with respect to A r and replace dS(X r — 
Xe + AX)/dX r with d5(X r — Xg + AA)/(i(AA). The differentiation d/d(AX) can be extracted to the outside of the integral 
because the integrand depends on AA only through this <5-function. 

As for (dQ / dt)[ lad \ with the aid of Eq.d53t. we arrive at the following expression: 



dg\ (rad) _ V 
dt/ t - T t 



dX' 



'nz(x')} / [e(x) (iy tr )d t + ly^d*) g^ ( x , z(x')) 



x=z{\)/ A 



2 Sl dX> 



(54) 



x'=z(\>) 
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Notice that the AA-differentiation in the second term operates only through z(X), not through z(X'). 



For (dQ/dt) 



(sym-S) 



dQ\ 
dt I 



we obtain a similar formula to the radiative part as 

(sym-S) 



dA'E[z(A')] 
dA' E(x') 



s(x)((v tr )a t 



d(AA) 



(V 0r HjG (syni - s) (^(A')) 
E[z(A)]G (sym - s) (z(A),x') 



(55) 



z'^zfA') 



where 



G {sym - S) (x,x') := G (sym) (x,x') - G {s) (x : x'). 
Using the following properties of G^ sym_,s ^ (cc, a;') (See Appendix AppendixlBlfor detailes.). 

G (sym - s) (2(A),z(A')) = G^ sym ~ s \z(X + A), z(A' + A)), 

" S) (x,x') = - 



_ G ( S ym-s) (a , )a/) 



° Msym-S), 



we find that the first term in Eq 

jq. (sym-S) 



vanishes to obtain 

q 2 d 



T t d(AA) 



dA'(s[z(A)]S[2(A')]G (sym - s) (z(A),z(A')) 



(56) 
(57) 

(58) 



(59) 



In Eq. d59) , the differentiation with respect to AA operates on both z(A) and z(X'), while, in the second term of Eq. (155) , it 
operates only on z(X). Since the part in the square brackets in Eq. ((59) is symmetric under the exchange of z(X) and z(A'), 
however, the contribution from the AA-differentiation through z(X') is the same as that through z(X). This is the reason why the 
factor 2 in Eq. d55) disappears in Eq. ((59). The detailed derivation is also presented in Appendix |B1 The result that the t- and <p- 
derivative parts in Eq. ((55) have no contribution can be expected from the fact that (dEjdt)t and (dL/dt)t have no contribution 
from the symmetric part of the Green function, which is shown by the balance argument of the corresponding conserved currents. 



We emphasize that, from the above expression, (dQ / dt) 



(sym-S) 



will not vanish in general. This means that the secular change 



of the Carter constant is likely to be also induced by the contribution from the symmetric part for the resonance orbit. 

It should be noted that we cannot divide the integration in Eq. (|59l into the symmetric and S-parts because both 
G( sym )[z(A),z(A')] and G^ [z(A), z(A')] diverge at the coincidence limit, z{X) — > z(X'). In order to treat them separately, 
we here introduce the point splitting regularization by displacing the orbits z(X) and z(A') as 



'(A) 



<(A) 



= *"(A) + 2*". 



«"(A')->^(A') 



*"(A') - 2^ 



with a small parameter Kl and the Killing field 



(Q^ cos C 



(60) 



(61) 



as the direction of displacement. ( is a parameter that specifies the choice of the Killing field, = T^/Tj and fl = T/T t . 
Under this point splitting regularization, we rewrite Eq. (l59l as 



dQ\ 
dt I 



(sym-S) 



lim 



>o T t d(AA) 

The potential functions, * (sym )(AA; e, C) and *( S )(AA; e, C), are defined by 



f (sym) (AA;f,() - * (s) (AA;e,C) 



* (X) (AA;e,C) 



+ 00 



dX' (e[z(A)]£[z(A')]G (x) [z+(A), z_(A')] 



(62) 



(63) 



where the superscript X is to be replaced to "sym" or "5". (Notice that, since S(x) is a function of r and 9, £[z±(A)] = £[z(A)].) 
It would be worth mentioning that the difference of the potential functions, >3> (sym_s ' :) (AA; e, C) := $ (sym) (AA;e, () - 
^(^(AA; e, £), is even in the sign flip of AA. To show this, we replace A and A' with —A and —A'. Under this replace- 
ment, the offset phase AA (and the Killing vector also flips the sign. Hence, we find lim e _>o <J/( sym-S, )(AA; e, £) = 
lim^o \J/( s y m_s ) (— AA; e, (). This symmetry guarantees (dQ / dt)[ sym S ^ = for AA = in general. This is consistent with 
the fact that the symmetric part does not contribute to (dQ/dt)t in the off-resonance case H, in which we can set AA = 0. 

Using the mode decomposed form of the radiative and symmetric Green functions given in Eqs. ((35) and ((37), the expres- 
sions in Eqs. ((54) and ((62) (with Eq. ((63)) can be recast into more practical expressions. We will show them in the following 
subsections. 
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B. Further reduction of (dQ/dt > f ad) : the radiative part 

Flanagan et al. lfl7ll have already written down a little more practical expression for the radiative part of the long-time averaged 
value of the rate of change of the Carter constant, (dQ / dt)[ lad \ in the gravitational case. Here, we briefly rederive the same 
result in the scalar case in a slightly different manner. 

Our first task is to rewrite Eq. (l54l in terms of the mode functions given in Eq. d2TT i. With the aid of the radiative Green 
function d35l l expressed as a sum of mode functions, we have 

( rad ) „2 r+oo +°° ' 



dt 1 1 167rzT t 7_; 

1 



ujlm 



^E E p 

m=-l 1 1 

d\'^[z(X')]} (-iu{V tr )+irn{V^) (nZf m [z(X)]) x 

UulmV U-oo J V d(AX) ) 

UPmU Klm? 1 f+OC d\'ntZ n W)}} (-ujQM+imQk.) - ^) <*f$T[*(A)]>; 



(64) 



where, for brevity, we introduced 7T^; TO (a;) := £(:r)7rj^ m (x), and also replaced the partial differentiations with respect to t and 
<fi with — ioj and im, respectively. We can rewrite 7r^ m [z(A)] as 

*ttM*)] = tUA-AA,A) e -^ T -™ T *> (65) 
with a function of two time variables, X r and Xg: 



(66) 

x=Az(X r ,X e ) 



where Az(X r , Xg) is the oscillating part of the orbit defined in Eq. dl9t . Since the function 3^; m (A r , Xg) is biperiodic with the 
frequencies T. r and Tg for X r and Xg, respectively, it can be expanded in the form of a Fourier series as 

oo oo 

ti m (X r ,Xe) = E E ~t nnT nMz- t(nrTrXrWeXe \ (67) 

72. 7 - — — OO fly — — OO 

where 

3Ln r nM ■= T- t dXr T~ r dXgt lm (X r ,Xe)e^™ M l 

A r Jo Ag J a 

By using Eq. d65l l with Eq. d67l ). we express the A'-integrals in Eq. ( l64l i as 

/oo OO OO ^ 

dX'tUMX')} = E E Y [ ^mn r n.He <T,rTrAA «(w - W ranrn( ) 

"°° TL r — — OO TLQ— — OO 

OO 

= E Z» mN (AX)S(u-u; mN ), (68) 



where we introduced 



and 



with 



N- 



w mnrnB := T t 1 (mT^ + n r T r + n e T ) = mfi^ + (j r n r + j ng)tt, (69) 
w m jv := + -Wfi, (70) 



^iroiv(AA) :— — 3/ m „ rnf ,(Wmn r n ( ,)e m '' TrAA , (71) 

(n r ,n e )6JV 



iV := {(n r ,n e )\j r n r + jgng = N} . (72) 
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Here it should be noted that N denotes an integer, and that we rearranged the order of summation over (n r , ng) in the second 
equality in Eq. d68l ). 

Owing to the delta function in Eq. d68l l. we can replace the Fourier integral in Eq. (164) with a discrete Fourier series. Hence, 



we can also replace uj in (jf i 



;(A)])a with uJ m N and then 



1 f T 

(A)])a = y 1 ™^^ j T dX H 3imn r n e ( w mJv)« 



i\(n r r r .+neTe-Nr)-in r T r AX 



E 3 b lmn r n a ( u mN)e 



i r A/2 

-m r T r AA_ / ^^ e «AT(n r >+ri6ije-A r ) 
A 7-A/2 



= ^Z** N (A\). 

Z7T 

In the second equality, we used the periodicity of the integrand with the period A. 
Finally, using Eqs. d68l l and (1731) . we can reduce Eq. (l64l to 



(73) 



dQ\ 
dt/ t 



(rad) 



E 



2ir 2 ^-r 1 |w m jv| 
1 



(0J m N(V tr ) ~ m(V^ r )) 



\ ZJ lmN 

(AA)| 2 



«?(AA)|< 



Zf-V(AA)tt(AA) 



Z^T(AA)yf°T*(AA) 



\ 



\ 



Z^(AA)y^(AA) 



UJ mNPmu; mN \ T uj mN lm\ 7 down/ A \ \ydown 



\v m NPmu mN \ \Pui mN lm[ 



32tt 2 Y l^mJvl 1 \au mN i m \ 2 

with := E£o £L=-i Ew=-oo and 

yU(AA) := ^^-^^(AA) = ^ ^ n^L^nX 

In the second equality we used the expressions of {dE/dt) t and (dL/dt) t (ll 



^m™(AA)yf r °™*(AA) 



,.T r AA 



(74) 



(75) 



d_(E 
dt 1 L 



6Att- 



E 



\UmN\ 



\ ZJ lmN 

(AA)| 2 , 



yaowi 

2 



(AA)| 



(76) 



The contribution from the radiative field in Eq. ((74) is nothing but the scalar analogue of {dQ / dt)^ 1 ^ obtained by Flanagan et 
al. 03 ■ Notice that the above expression for the radiative part is also valid for the off -resonance case once we interpret that the 
summation over N means the summation over independent frequencies in this case. 



C. Practical formula of {dQ/dt)\ Byjn ': the regularized symmetric part 

To obtain a more practical expression for (dQ / dt)^ ym ~ S ' , we rewrite the potential functions of the symmetric and S'-parts 
given in Eq. (lo*3l , respectively. The symmetric part can be treated in a similar manner to the radiative part discussed in Sec. lIIIBl 
though we cannot completely separate the averages over A and A' at the level of each mode. Substituting Eq. ( |37| l (with Eq. ( l38| >) 
into the potential function of the symmetric part (|63) , we obtain 

i r°° 

/oo 
dX' (^e-^-^^-^g^iAziW- AA), Az(A', A' - AA)]) , (77) 
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where 

Since ff^im^ [Af (A, A — AA), Az(A', A' — AA)] is a periodic function with period A, we find that the A'-integral in Eq. (T7Tb 
produces 5(uj — uj m N), and therefore we can discretized the Fourier integral as 

*^ m >(AA; e,0 = ^ e^^^i ^dA' ^e~* JVT ^~ A '^ ( ^™] TO [Af (A, A - AA), Az(A', A' - AA)]^ , (78) 

where, e\ = ecos£, and £2 = esinf. In the above equation, the long-time average over A can be replaced by an average over 
one period because the part in the angle brackets is periodic with period A. For later convenience in subtracting the the potential 
function of the S-part, #( s ) (AA; e, (), from ^( s y m ) ( AA; e, C), we rewrite Eq. d78]l by rearranging the order of the summation 
with respect to I and to as 



*^ m '(AA;e,C) = £ J2 e- in ^ N+ ^ m ^^ (AX) (79) 

N= — oo m— — 00 



with 



1 00 1 /-A pA 

^ m) (AA) :=^r E a? / dX J ^'e^^'^LlMA, A - AA),^', A' - AA)]. 

l—\m\ ^ ^ 



(80) 



We interpret Eqs. (l79l and (l80t as the Fourier expansion in the (ei , £2 )-space and the Fourier coefficients, respectively. The point 
is that each (to, N) Fourier mode is finite, which can be calculated numerically. 

Next, we derive "f^* 1 (AA; e, £). Using half the squared geodesic distance between x and x', a(x,x r ), the S'-part Green 
function defined in the local convex neighborhood of the particle's position is written as 

G (s) (x,x') = -— [U(x,x')S(a(x,x'))-V(x,x')9(a(x,x'))} (81) 

with 

U(x,x') = l + 0(\x-x'\ 3 /L 3 ), (82) 

where L is the scale of curvature of the background spacetime, and V(x, a;') is a regular function in the coincidence limit x —> x' . 
(See, e.g., Ref. yfl for more concrete definitions of these quantities.) The integrand containing V(x, x') contributes as the finite 
value to ^^(AA; e, C) since A' integral has its support only in a short interval around A' w A. In the coincidence limit, this 
interval shrinks to zero, and thus this contribution also vanishes. Then, we need to consider only the term including U(x, x'). 
The argument of the <5-function is expanded, abbreviating the argument AA, as 

a[z + (X), Z-(X')} = \ gilv [z{\)\ K(A) - z^(X') + eC} K(A) - z"(A') + ef} + 0(e 4 ) 

= \a^[z(X)] {£[z(A)K(A)<SA + ee 1 } {E[z(X)]u"(X)5X + e?} + 0(e 4 ), (83) 
where we introduced new variables 

<5A:=(A-A'), A:=i(A + A'), (84) 

and assumed SX = 0(e) since the solution of a = for SX is so. We stress that 0(e 3 ) term is absent since a is unchanged under 
the transformation e — > — e and <5A — > — SX. With the aid of the above expansion for a, expanding both £[z(A)] and £[z(A')] in 
terms of e as well, we eventually obtain 

* (s) (AA;e,C) = + 0(e) (85) 
e 

with 



4 ^ A Jo x/(g^ + u^)^e 



(86) 

x— z(A) 
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where we replaced the long-time average over A by the average over one period because the integrand is periodic with period A. 
Note that i/j(Q depends on £ through ^ . This is the only S-part contribution that remains in the limit e — > 0. Since (g^ + u il u u ) 
is a spacelike projection operator, the expression inside the square root is positive semi-definite. 

To subtract the S'-part of the potential function from the one of symmetric part mode by mode, we calculate the (to, 7V)-modes 
of ^(S)(AA;e,C) as 



4tt 2 



TI /'TI 

dei 



~4tt 2 



n 

2^2 



d eae <n(a^+«m)«r(S)(AA;e,C) 



17 dd de 2 e 8 °( £lJV+e2m ) 



cosC 



deie 4ei0(Ar+mtanC V(C) 



(-1)' 



N cos £ + to sin £ 



sinC 



de 2 e "2f2(WcotC+m)^^ 



sin(7TTO tan £) VKC) 
iV cos C + to sin £ 



(87) 



This expression is a double integral with respect to A and £, and it is finite. Since ip(0 is independent of to and TAT, once we 
obtain tp(() as a function of ( by performing the A-integral, the remaining integral over ( for each pair of to and TV is just a 
single integral. 

We stress that subtracting Eq. ((87) from Eq. ( |80l is a variant of the so-called mode sum regularization 01 developed in the 
context of self-forces, but our proposal is suitable for the decomposition in terms of the spheroidal harmonics. We do not have 
to reexpand the perturbation using the spherical harmonics, which is usually utilized. After subtracting this S'-part contribution 
from Vl/^J" (AA) mode by mode, one should be able to take the summations over to and A^ without any divergence. The 
convergence of these summations will be even accelerated by fitting the asymptotic form as usual. We will report the result of 



explicit convergence test of ( AA) 



^> mN {l±.X) in our succeeding publication. 



IV. SUMMARY AND FUTURE DIRECTIONS 



In this paper, we have derived the long-time averaged rate of change of the Carter constant, (dQ/dt) t , for a self-interacting 
scalar charge within the adiabatic regime, which is also applicable to a resonant inspiral. The essential point is that, by contrast 
to off -resonance orbits, a resonance orbit depends on the offset phase, AA, which is the difference of the time when the r- and 8- 
oscillations reach their minimum values. While this offset phase is identically zero for a off-resonance orbit, it takes a non-zero 
value for a resonance orbit, and also evolves slowly compared to the time scale of the orbital motion. Therefore, we need to 
care about the evolution of not only the three constants of motion but also AA. The existence of a non-zero AA requires a 
direct computation of the self-forces acting on the charged particle even for calculating the averaged rate of change of the Carter 
constant. This means that (dQ/dt) t for a resonance orbit is no longer determined solely by the radiative field $( rad ) (x) alone, 
and that we need to consider the contribution from the symmetric field, $( s y m ) (x) as well. 

As for the radiative part (dQ / dt)\ iad \ we derived a practical expression in a similar manner to an off-resonance case with 
slight modification, and found that our result is equivalent to the scalar analogue of the long-time averaged rate for the gravi- 
tational case shown in Ref. Il7ll . On the other hand, the symmetric part is more complicated than the radiative one because of 
the singular behaviour of the symmetric field at the location of the particle. We proposed a method of subtracting the singular 
behaviour, so-called S-part, from the symmetric part via the point splitting regularization in the Killing directions. The key 
point is that it is compatible with the Teukolsky formalism, which makes the calculation much simpler than the direct self-force 
calculation. While there is an expectation that (dQ / dt)[ sym S ^ might eventually turn to be supressed, which is suggested by 
post-Newtonian analysis [28], this conjecture should be tested in the fully general relativistic analysis. Since the numerical code 
that calculates the long-time average of the constants of motion for a generic off -resonance orbit around the Kerr black hole is 
well developed 129113011 . there is in principle no obstacle to numerically test the conjecture. The code is now under development 
and the result will be reported in our succeeding publication. 

To avoid unnecessary confusion over various works, we comment on the following two points related to our current works. 
The first point is the difference between the long-time average and the phase space average for a resonance orbit, discussed in 
Refs. |0, [Ot]. Different from the long-time average, the phase space average for a resonance orbit means that the averaging is 
further performed with respect to the slowly evolving variable AA that is an important control variable of the resoance orbit. 
The phase space average eliminates all dependence of A A both in (dQ/dt)[ lad ^ as well as (dQ / dt)[ bym , thus it is likely to 
miss the important contribution to the accumulated orbital phase correction that would exceed 0(1) radians. (See discussion in 
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Appendix [A]) For this reason, we believe that only the long-time averaged rate of the change of constants of motion can be a 
good candidate to obtain the adiabatic evolution of a resonant inspiral with phase accuracy better than O(l). 

The other poinst is that our new method does not reduce the importance of developing the direct self-force calculation that is 
recently under a rapid progress. We think that these two methods are complimentary. Relatively high precision will be required 
in evaluating the effect that may cause an orbital phase error enhanced by an inverse power of the mass ratio. Numerical 
evaluation of such quantities might be a tough problem if we cannot take advantage of the Teukolsky formalism that reduces the 
problem into the one to solve one-dimensional radial functions. Even if achieving high accuracy numerically is not so difficult, 
our method will give reference numbers which are useful to check the accuracy of the results obtained by direct self-force 
calculations. 

Before closing this paper, we would like to mention how to extend our method to the case of the gravitational radiation 
reaction. Since our method is based on the mode sum regularization in which each mode contribution is gauge invariant, we 
do not have to care about the gauge adjustment between the symmetric part and the S'-part. Even for the gravitational case, 
the contribution from the S'-part evaluated in the Lorenz gauge is almost identical to the scalar case discussed in this paper. To 
calculate the contribution from the symmetric part, we can use the standard Teukolsky formalism with the aid of the established 
method of reconstructing the metric perturbation from the master variables (HU-HH]. The metric perturbation reconstructed in 
this manner is given in the so-called radiation gauge, in which a singularity runs from the particle's trajectory in the radial null 
direction. However, almost all regularized trajectories shifted in the Killing directions do not cross this singularity. Therefore 
we expect that the presence of this radial singularity will not cause any trouble. Although we need further study for practical 
implementation, we believe that our method can be extended to the gravitational case as briefly described above, 
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Appendix A: The possible pattern for the long-term orbital evolution crossing the resonance 

The final goal of this project is understanding the long-time orbital evolution crossing the resonance in the gravitational case. 
Toward this ambitious goal, this appendix is dedicated to quantitatively discuss how we can understand the long-term evolution 
of the resonant inspiral, as a future application of our formalism. In this appendix, we consider the gravitational case: a point 
particle with the rest mass fi moves along the quasi-resonant orbits in the Kerr black hole with the rest mass M and the azimuthal 
angular momentum aM. To apply the following discussion to the scalar model, we just need to replace fi to (q 2 /fi) in the results, 
where q is the scalar charge of apoint scalar particle. 

As was discussed in Ref. SIM], as far as we consider the first order self-force in the osculating orbit approximation, the long- 
term orbital evolution is basically determined by tracing the evolution of slowly changing variables, i.e., the constants of motion 
P, and the orbital frequencies as their functions, T a (P). We only need to know the non-oscillating part of the rates of change 
of P and T a (P) for a generic orbit as the corrections due to the self-force. In particular, we only need to know the former if it 
is allowed to neglect O(l) error in the orbital phase. This observation is essentially equivalent to the more systematic two-time 
scale analysis given in Ref. @]. 

Near the resonance, however, the adiabatic evolution of the orbit cannot be described by tracing the constants of motion P 
because the self-force also depends on the offset phase AA, the relative phase between r- and (9-oscillations, as we have discussed 
in the main text. Neglecting the part that oscillates in the orbital time scale, the time derivatives of P are determined as functions 
of P and AA, which we denote by P (P , AA) . At the same time, however, A A also evolves unless the exact resonance condition 
is satisfied, and the resonance orbit effectively comes back to the off-resonance orbit when the offset phase starts to change 
rapidly. We examine possible patterns of evolution of {P, AA} in this section. 

Suppose that the resonance condition is only approximately satisfied, but its deviation is small. In this situation, P may still 
depend on A A. Moreover, we are also allowed to replace P in P(P,AX) by the fixed value at the exact resonance point, 
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which we denote I* ea . This can be understood in the following manner. The typical time scale for crossing the resonance point 
is estimated as r ros ~ 0(CIaWM //i). During the resonance crossing, the net change of the constants of motion from I* es is 
estimated as AP ~ I l T rcs ~ 0(fl^fi/M) 0,0]. This means that the correction to I* es is always suppressed by the factor of 

0(y/fI/M) in resonance crossing. Based on this observation, we can estimate the error in P that associates with the replacement 
of its argument from P to 7* es as 

j f (log AA)) = ± (log4 s (/\AA)) +0 (n* (£:) 3/2 ) ■ (Al) 

After passing the resonance, this error term causes the shifts of the constants of motion, which will be evaluated as A/| rr /7 l ~ 
O(fifM). However, the net phase error due to these shifts, which is simply given by A/| rr times r ra d is at most 0(1). This 
phase error is small compared with what we are interested in here. Thus, we conclude that the evolution of the constants of 
motion during the resonance crossing is given by P(I* es , AA). 

Next, we discuss the evolution of the offset phase, AA. In the situation considered here, there is no longer exact common 
resonance period A. Instead, we introduce the approximate common period A by 

A:=j e A e = ^. (A2) 

During the period A, the offset phase is shifted by 

<5(AA) = A - j r K = %^AT, (A3) 

Z7T 

where AT :~ jeT r — j r Yg. When AT — 0, we have the exact resonance. In the adiabatic regime, we are also allowed to 
replace the time derivative of AA to the averaged value over one period A. As a result, the evolution of AA is described as 

/ dAX\ _ X 5(AA) _ _AT_ 

\ dt /- Tt A " j e r t r r - (A4) 

Now, we investigate the evolution of AA around a resonance point in more detail. For this purpose, focusing on the time 
derivative of (dAX/dt) t , we consider the two cases, separately, depending on whether or not there exists AA C such that 

(See also Ref. l35ll for the classification of resonance orbits.) The latter case has been already discussed in Ref. HI] and is 
essentially the same as the transient resonance analysed in Refs. lfl3i[l7ll . In this case, AA varies of 0(1) on the time scale of 
0(Qj, I y/Ji) around AT = 0, during which P evolve by 0{ s /Ji) relatively. The effect of the resonance alters the evolution of P 
which leads to the phase corrections of 0{- s fJXl < j,AT), where AT is the time interval between the resonance and the plunge. As 
typically AT scales like (oc l//i), the phase corrections due to this effect can be significantly large (oc 1/ -JJi) if the inspiral in 
the extreme mass ratio regime, M /i. 

A more interesting situation may happen if AA C exists. This case corresponds to the so-called sustained resonance lfl3tl . When 
AA C is sufficiently close to AA C , we obtain, with Taylor expansion around AA C , 

j t (^} = -^ 2 (/ 4 )(AA - AA C (!')), (A6) 

where 

Since the time scale l/\vo\ scales like (oc 1/y/Ji), it is likely to be much shorter than the radiation reaction time scale, which 
scales like (oc l//i). If w 2 > 0, the above system behaves like a harmonic oscillator with its potential minimum varying 
adiabatically. In this case, AA oscillates around AA C with the amplitude oc l/y/w. This evolution may last as long as the 
condition w 2 > and the existence of AA C are maintained. Therefore, in principle, the resonance may last for the radiation 
reaction time scale. This time scale is much longer than we naively expect when the possible presence of AA C is neglected. If 
vd 2 < 0, the value of AA close to AA C is disfavored. The time scale for the variation of AA is again oc 1/ y/j2 as in the case of the 
absence of AA C . Hence, the corrections to the phase evolution are also of the same order as before. Although the post-Newtonian 
analysis lfl3ll and the gravitational wave fluxes from a resonance orbit [fl7,j suggest that the sustained resonance may be absent in 
inspirals around the Kerr black hole, or exist only in very restricted orbital parameter regions even if it is achieved, the existence 
of the sustained resonance in inspirals into the Kerr black hole is still an open question. 
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Appendix B: Derivation of Eq. l |59l 

First, to show Eqs. (l57l i and (l58l . we clarify the properties of z M (A) and G^ sym ~ s ' (x, x'), respectively. For resonance orbits, 
r(A) and 6(X) are periodic with period A(= j r A r = jgAg) as shown in Eq. (|7}, while t(X) and <fi(X) have secular evolution 
terms proportional to A in addition to oscillatory parts with period A as Eqs. $1% and (1141 1. Taking account of these facts, we 
find that, under the transformation A — > A + A, z M (A) transforms as 

z^(X + A) = z"(A) + [Tttfa + T^J A = *"(A) + T t A^(0), (Bl) 

where £ M (C) is defined in Eq. (loTt . This means that the trajectory is displaced in a Killing direction. 

As for G^ sym_S ' ) (x, x'), we find two properties. Due to the stationarity and axisymmetry of the Kerr spacetime, the Green 
functions are invariant under the displacement along an arbitrary Killing vector field: 

G(x + M(Q,J + M(0) = G(x, x'), (B2) 

where A is an arbitrary constant. From Eqs. (IBU and (1B2K we can derive Eq. ( fS7b as 

G(z(A + A),z(A' + A)) = G(z(A) + T t A£(0),z(A') + T t A£(0)) 
= G(z(X),z(X')). 

The relation dB2b simply tells that the Green functions depend on t and t' through t — t', and on and <fi' through <f> — <f)' . Based 
on this understanding, we find that Eq. d58l holds. 

Next, we show that the first term on the right hand side of Eq. (l55l l vanishes by using Eqs. d56*l l. d57b , and 08] ). We focus only 
on the t-derivative part, but one can show that the ^-derivative part vanishes in the same way. By using the property d57l i. we 
can replace the long-time average over A to the average over one period A as 



dX'nz(X')](jX(x)d t G(x,z(X'))) x=zWj 

i r T r°° 

= J T dX J jxn<mi4nidtG(x,x>)} x=z{x]x , =z{y) 



l-T 

N-l 



lim 

jv->oo 2NA 

n=-N 



-i /-(n-rxjn poo 

— Y, / d\ dA'E[z(A)]E[z(A')][a t G(x,x')] x=z(A) ^ =z(v 

N-l »A poo 

n^oo 2NA E J o ^/_ oo ^'S[z(A)]S[^(A')][9 t G(x,x')] x= , (A)iK , = ^ ) 



pA poo 

A 



i / dA / dA'Sl^^lS^CAOl^G^x')]^),,,^). (B3) 

'0 J —oo 



In the third equality, we transformed (A, A') to (A — nA, A' — ?tA), and used Eq. $5% and E[z(A — nA)] = E[z(A)] (note that 
£(x) is a function of r and 0). We also rewrite the integrals over A and A' as 



j f dX I d\'E[z(\)]E[z(\')] [d t G(x, x')] x=zWtXl=z{xr) 



1 /-A /-(n+l)A 

t dX «iA'E[z(A)]SKA')][9 t G(^')] I= , ( x),^,(v) 

1 °° /■ — (n-l)A <-A 
' v „ J -nA JO 



A 

n— — oo 



-nA 

rfA / A dA'E[z(A)]E[z(A')] z')L=.(A),*'=*(A') ■ (B4) 



Further, for the Green function G < - sym s ' (x, x'), by exchanging the labels in the last line of Eq. ( IB4I ) as A O A' and x O x', we 
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obtain 



A 



dA / dA'E[z(A)]E[«(A')] d t G ( - sym ~ s \x, x') 

-oo JO L 
1 

A 



2(A), x'= ^(A') 

dA / dA'£[z(A)]E[.z(A')] ^G (sym ~ s) (x',x) 







x=z(X) 1 x r —z(X') 



A 



dA / dA'E[z(A)]E[z(A')] -d t G( sym - s )(^') 
J —00 

where we used Eqs. (156l l and (158l l in the last equality. Therefore we find 



dA'E[z(A')] ( h(x)d t G {sym - s \x,z{\')) 



x—z(X),x f —z(X f ) 



= 0. 



x=z{\)/ x 



(B5) 



(B6) 



As for the AA-derivative part, we transform 

d 



£(*') 

A 





dA' 

OO 

= i / 



x S[z(A)]G^ m - s )(z(A) ja; ') /; 



z'=z(A') 



E(x) 5 ^r{E[z(A / )]G^-«(x,«(A'))} 



(B7) 



x— 2(A) 



where we replaced the long-time average with respect to A with the average over one orbital period A as shown in Eq. ( lB3l l. and 
exchanged the labels as A H A' and x <-> .t', and the domains of integration with respect to A and A' in the same way as in 
Eq. dHl). Using Eq. dB7]), we find 



dA' 



E(a: ' ) d(AA) {n^)\G (sym - S) (^)^') /; 



x' — z(X') 



1 

2A 



1 

2A7o 
1 d 



dX / dA' 



S(x ' } d(lA) { S[z(A)]G(Sym " S)(z(A) " T ' } } 

£(*)^^{eKa')]g^- s W(a'))} 



,(A') 



X— 2(A) 



dA / dA'- 



E[z(A)]E[z(A')]G^ s )(z(A),z(A')) 



2 d(AA) 



d(AA) 

dA' (s[z(A)]E[z(A')]G (sym - s) (2(A), z(A')) 



(B8) 



Removing the t- and 0-derivative parts from Eq. (I55t and using Eq. (IB8K we finally obtain Eq 
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